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Abstract 

Form factors are calculated in the point form of relativistic quantum mechanics 
for the lowest energy states of a system made of two scalar particles interacting via 
the exchange of a massless boson. They are compared to the exact results obtained 
by using solutions of the Bethe-Salpeter equation which are well known in this case 
(Wick-Cutkosky model). Deficiencies of the point-form approach together with the 
single-particle current are emphasised. They point to the contribution of two-body 
currents which are required in any case to fulfil current conservation. 
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1 Introduction 



Calculations of form factors often retain in a first approximation the single-particle current 
(impulse approximation). In most approaches, this contribution has to be completed by 
at least two-body currents, especially to fulfil current conservation. This also holds for 
relativistic approaches where a covariant calculation of the single-particle current does not 
necessarily imply that physics is properly accounted for. For each relativistic approach, it 
is therefore important to test the degree of validity by comparing its predictions to a case 
where an exact calculation can be performed. Neglecting vertex and mass corrections, as 
usually done on the basis that they are partly incorporated in the physical inputs, such 
an exact calculation is generally believed to be provided by the Bethe-Salpeter equation 
with an appropriate interaction kernel. 

A particular case of interest is the Wick-Cutkosky model 0, where the Bethe-Salpeter 
equation is solved in the ladder approximation while the interaction kernel results from 
the exchange of a scalar zero-mass boson between two distinguishable scalar particles. 
Solutions can be obtained relatively easily due to an extra hidden symmetry. Using 
the expressions of the Bethe-Salpeter amplitudes, calculations of form factors can be 
performed for the lowest bound states that we intent to consider here. Contrary to other 
approaches, there is no need to add two-body currents. The contribution of the single- 
particle current is sufficient in the present model to ensure current conservation, which 
can be checked with the expressions of the matrix elements of the current]^. This model 
was used by Karmanov and Smirnov as a test of the description of form factors in the 
light-front approach for systems composed of scalar particles with small binding energy 
1^. They used the non-relativistic expression, — for the binding energy, which differs 
from the exact one, but this does not seem to affect their conclusion. Interestingly enough, 
the comparison of the exact calculation with the non-relativistic one does not show much 
difference up to ~ 100 where m is the constituent mass. Beyond, relativistic 
corrections with a log character slowly begin to show up. 

We propose to make a similar test for the point form of relativistic quantum mechanics, 
which is one of the forms proposed by Dirac, beside the instant and the light-front forms 
0. This form, which is much less known than the other two, has been developed in 0] 
and recently used for a calculation of form factors of the deuteron and the nucleon ^ . 
In both cases, the form factors decrease faster with than the non-relativistic ones. In 
the first case, the discrepancy with experiment tends to increase while, in the other one, 
it almost vanishes. We will not extend in this letter on the questions that are raised by 
this last observation. We will only present a few results which, by themselves, are quite 
significant. A more complete analysis will be presented elsewhere fl^ . 



In the present study, we calculate the ground-state (/ = 0) form factor as well as a 
transition form factor to the first radial excited state. Different couplings are considered, 
corresponding to states weakly, moderately and strongly bound [a = 1.0, 3.0 and 27i 
in terms of the QED coupling). The last value implies a total zero mass of the system, 
which is an extreme (unphysical) case, nevertheless interesting to look at, too. Momentum 

^This result together with some of the expressions used in the present paper for the Wick-Cutkosky 
model are part of a work in preparation. A preliminary presentation was made in B . 
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transfers, Q^, up to 10 times the constituent mass squared, m^, are considered. These 
different cases will provide a sample of results which are significant enough to test the 
validity of the single-current approximation in the point-form approach and give insight 
on the results presented in refs. |^. 



2 Expression of the single-particle matrix element in 
different approaches 




Figure 1: Representation of a virtual photon absorption on a two-body system with the 
kinematical definitions. 

The contribution which we are interested in is shown in Fig. 1. The general expression 
of the corresponding matrix element between two states with / = 0, possibly different, is 
given by 

^/2Ef 2E, (/I J'^K) = (Pf + Pt) + F,{q') (1) 

where = —Pf^. Current conservation imposes constraints on the form factors Fi{q^) 
and F2{q^). For an elastic process, F2{q^) has to vanish but this result automatically 
stems from symmetry arguments alone. It does not imply that current conservation holds 
at the operator level, as it should. For an inelastic process, the following relationship has 
to be fulfilled: 

Fi(g2)(Mj-M2) + F2(g2)g2 = 0. (2) 
Form factors using the Bethe-Salpeter amplitudes. 

For the model under consideration here, the general (and exact) expression of the matrix 
element of the current, which reduces to a single-particle one in this case, can be written 
in terms of the Bethe-Salpeter amplitudes, Xp{p)i 

^2Ef2E, (/|J^|2) = 

{Pf+Pt-2p^) {p'-m^)x,^{lp.-p). (3) 

For the Wick-Cutkosky model, the Bethe-Salpeter amplitudes take the form of a relatively 
simple integral representation, i.e. for the lowest energy state with a given orbital angular 
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momentum /: 

"1 ^?„(z) YriP) \p\ 



1 [m? — jP"^ — — z P ■ p — ie)"+2 ' 

with n = I + 1. In this expression, gn{z) is the solution of a second order differential 
equation ||l|, ||, that can be solved easily. 

Form factors in the non-relativistic limit. 

Opposite to the full relativistic calculation, a non-relativistic one can be performed. Using 
wave functions that are solutions of a Schrodinger equation, general expressions for both 
elastic and inelastic form factors can be obtained: 

It can be checked that the above form factors verify the current conservation condition, 
Eq. (0), provided the interaction is local. The second form factor vanishes in the elastic 
case. 

Form factors in the point-form approach. 

It has been shown |^ that a calculation of form factors in the point form approach could 
be performed relatively easily by using standard wave functions obtained from a mass 
operator of the form 

M = Mfree + M.^f (6) 

This includes a large class of wave functions, since the sum of the kinetic and potential 
energies appearing in the standard Schrodinger equation can be identified with the op- 
erator (up to a factor). This holds for a two-body system and provided the energy 
E is appropriately redefined. Therefore, wave functions entering the non-relativistic ex- 
pressions of Eq. d^) can be used. How the matrix element of the single-particle current is 
calculated has been described in ref. 0. Instead of using an expression where appropriate 
boosts have to be performed, we here give one whose Lorentz-covariance is explicit: 



2Ef2Ei (/|J^|z) = J2M/2M, 



(27r)3 

X J d'^pd'^pf d'^pi drjf drji 6{p'^ — rn?) 5{p'j — rn?) 5{p1 — m^) 
xe{Xf pf) e{Xf ■ p) e{\ ■ p) e{\ ■ p,) 6\pf +p- Xfr]f) 6\pi +p- A,r/,) 

Hi^r) ^{pf+p?{p^+pf (p'f+pf)- (7) 

In this expression, all quantities are Lorentz-invariant ones (except obviously for the 
current which behaves as a 4- vector). The auxiliary variables, rji and rif, have been 
introduced to make the covariance manifest. When these variables are integrated over, 
they give rise to 3-dimensional 5(...)-functions 

S{p^ + p-^ip■+p'')) and 5(p/+p-^(p;+/; 
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that are essential relations in the point-form approach. The Af and Aj 4-vectors are 
unit vectors proportional to the 4-momenta of the total system in the initial and final 

states, Af = and Aj = jjj. They can be expressed in terms of the corresponding 

velocitie^, A° = (Vl — v^)^^ and A = v{\/l — v^)~^ (c=l)- In the cm., it can be 
checked that the wave function only depends on the relative momentum of the 

two particles. On the other hand, it can also be verified, by direct integration or after 
performing a change of variable, that the current of the system under consideration is 
given by ( , (J)) = i^yv), in agreement with the standard normalisation of the 
wave function 



(27r)3 

Expressions of the form factors, Filq"^) and -F2(g^), can be obtained in any frame. With 
an appropriate change of variables, they can always be expressed in the forms they take 
in the Breit frame (here defined hj v = Vf = —Vi), where they may be simpler. Using 
auxiliary quantities, Fi{q^) and ^2(5^), they read: 



F,{q') J2Mf 2Mi = Fi{q') {Mf + M,) - F2{q') {Mf - M, 



with 



In the above equations, the velocity defined in the Breit frame, is related to the mo- 
mentum transfer by the relation v"^ = Q2^(j[.j^_|_j|^')2 • The (Lorentz-) transformed momenta 
are defined as: {p^ ,p^ ,p^)tf = {p^,P^, y^^) ^ind {p^ ,p'y ,p^)ti = {p^,p^, ^7^^), together 



with Cp = y/ m'^ + p'^ 



Wave functions and analytical results. 

For the wave functions of the ground and first excited states, (f){p) and (f)*{p), we use 
solutions obtained with a Coulomb-like potential. 



= \/47r- 



(k^ -|- p2)2 ' 

<^*(P^) = v^(-^^7-i)l(-"-p")' (10) 

where = — |M^, the total mass M being that one obtained from the Bethe- 
Salpeter equation for the ground state. It has been shown [0 that the spectrum of 
the normal states for the Wick-Cutkosky model, which has the same degeneracy as the 
Coulomb potential, could be reproduced with a 3-dimensional equation and an effective 



^The quantity, v, here refers to the usual velocity 3-vector, v, a notation that differs from the one 
employed in ref. where it represents a 4- vector which corresponds to our 4- vector A** . 
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interaction of the Coulomb-typeQ Due to the appearance of in such an equation, 
the energy redefinition mentioned previously is not even needed in the present case. The 
wave functions of Eqs. ([T0|) should therefore be a good zeroth order approximation for 
our study, including the extreme case = 0. Accordingly, we assume k*^ = 0.25 = 
0.25 {maeff/2y, which differs from the Bethe-Salpeter result, k*^ = — ^M*^, by a few 
percent. With these wave functions, the form factors can be calculated analytically. We 
checked that the main features evidenced by our results were insensitive to this choice 
by using a different (numerical) wave function [T^, more in the spirit of the point-form 
approach, i.e. obtained directly from the linear mass operator, Eq. (^. 

In the non-relativistic case, the elastic and inelastic form factors are respectively given 
by: 

^-(-^^)- (.^ + gW ^'^'^^^'^ 

Taking into account that Mj — Mf = one can verify that the condition of current 
conservation, Eq. (H), is fulfilled. 

In the point-form approach, the elastic form factor, written in a way that resembles the 
non-relativistic one, reads: 

Fi{q' = -Q') = - } ,,.4 . F,{q' = -Q') = 0. (12) 
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Interestingly, the factor 1/(1 + -^p) which multiplies the quantity in the denominator 
of Fi{q^) is the same as the one sometimes introduced by hand in order to account for the 
Lorentz-contraction effect (see discussion in ref. Contradicting asymptotic results 

(in QCD for instance), the range of validity of this recipe is limited to small Q^. Curiously, 
part of the extra factors also look like a Lorentz-contraction effect. 

As for the inelastic form factors for a transition from the ground- to the first radially 
excited state, they are most easily expressed in terms of the quantities, Fi{q^) and F2{q^), 
which can also be calculated analytically: 

Fi g2 = _ 72 \ ^ J' 1 + v^) 1 - v^, 

(9ft"^ + (16 m"^ — IOk"') + v^K^y 

= (0 2^ 2nyr?n"l^ 4 2^3 (1 + (1 - (13) 
(9^"^ + (16 — IOk"') + v^n'^y 

where is defined after Eq. (^. These expressions generalise those for the non-relativistic 
case, Eq. (^) (it is reminded that in this limit = 16f^m^). In contrast however, they 
do not allow us to fulfil current conservation, Eq. (Q). 

■^The spectrum of normal states of the Wick-Cutkosky model is reproduced within a factor 2 by an 
equation of the type [4(to^ +P^) — A/^] (/)(r) ~ 4m </'(?'), even for cte// oo. 
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3 Results 



Table 1: Elastic form factor, Fi(g^), for the ground-state. Units for E and k are the con- 
stituent mass, m. Results are presented for different couplings and, for each of them, for 
different approaches: Bethe-Salpeter equation {B.S.), non-relativistic calculation (N.R.) 
and point-form approach {P.F.). The value of a referred to in the table corresponds to 
the Bethe-Salpeter equation while the coupling for the "non-relativistic" model, cueff, is 
chosen in such a way to reproduce the binding energy for the ground state. 



Q /m 


O.Oi 


n 1 
O.i 


1 n 
i.O 


inn 
iO.O 


a = 1 










E = 0.0842, = 0.0824 










B.S. 


0.984 


0.856 


0.309 


0.137-01 


i V .it. 


u.yoo 






u.ioo-ui 


P.F. 


0.984 


0.853 


0.298 


0.097-01 


a — 3 










E = 0.432, ^2 = 0.385 










B.S. 


0.996 


0.962 


0.705 


0.139 


N.R. 


0.996 


0.968 


0.740 


0.145 


P.F. 


0.994 


0.948 


0.620 


0.056 


q; = 27r 










E^2.0, ^2 = 1.0 










B.S. 


0.998 


0.983 


0.848 


0.339 


N.R. 


0.997 


0.987 


0.884 


0.378 


P.F. {E = 1.90) 


0.613 


0.397-01 


0.110-03 


0.124-06 


P.F. {E = 1.95) 


0.181 


0.142-02 


0.191-05 


0.196-08 



In Table 1, results are presented for elastic form factors corresponding to small {a — 1), 
moderate {a — 3) and strong binding {a — 27r). In the last case, the elastic form 
factors in the point form vanish at finite and results are actually those obtained when 
approaching the limit a 27r, with E = 1.90 m and E = 1.95 m. For these energy values, 
the non-relativistic results, to which the previous ones may be compared, are essentially 
the same as for E = 2.0 m, given in the table. 

One immediately notices that the non-relativistic results are very close to the exact 
ones for all cases, including the extreme case where the total mass of the system is zero. 
This agreement extends up to Q"^ — 100 with an error of 20% for a = 1 and 50% for 
a — 27:. The discrepancy at small momentum transfers is typically of the order of y^^- 
Most probably, it can be traced back to the wave function used in the non-relativistic 
calculation or to the electromagnetic single-particle current. These ingredients do not 
fully account for corrections due to factors ^ (in the potential for instance) or for the 
field-theory character of the Wick-Cutkosky model. A major lesson of these results is 
that relativistic effects are not necessarily important and that most probably, for any 
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non-relativistic calculation of form factors, there exists a covariant calculation (the exact 
one in the present case) which gives very close results over a large range of momentum 
transfers. However, this does not mean that this calculation involves all relativistic effects 
and is physically relevant, especially with respect to current conservation. An example is 
provided by the deuteron electro-disintegration near threshold in the light-front approach. 
Two covariant calculations of the transition form factors have been shown to be very 
close to the non-relativistic ones up to = 10 (GeV/c)^ [0, completely missing the 
contribution due to the pair term whose relativistic character is well known. In the non- 
relativistic approach, this contribution, which is required to fulfil current conservation, 
is essential to account for experiment in the low momentum transfer range. How this 



contribution appears in the covariant light-front formalism was shown later on [|I4 



The comparison with the point-form results evidences a discrepancy that increases with 
the momentum transfer as well as with the coupling strength. It becomes especially large 
when approaching the extreme case where M = 0. Two features are worthwhile to be 
mentioned, that stem from examining the analytic expressions of the form factors, Eq. 
(p!^). First, there is a contribution to the squared-charge radius which varies like j^, as 



it was found numerically in ref. p]. Second, the form factor drops more quickly with 
than in the exact or in the non-relativistic calculations, roughly like ^ instead of 
This can be seen in Table 1 for E = 1.90 m and E = 1.95 m and is in agreement 
with what the examination of the Born amplitude reveals. The effect becomes especially 
sizeable when approaching = 4 and beyond. 

Table 2: Inelastic form factors, Fi{q^) and F2{q'^), for a / = — / = 0* transition, 
and a = 3,Ei = 0A322m,K^ = 0.385 m^E/ = 0.1036 m, = 0.101 m^ for the exact 
calculation and 0.096 m^ for the non-relativistic one. 



QVm2 


0.01 


0.1 


1.0 


10.0 


B.S. 












-0.032-01 


-0.298-01 


-0.145-00 


-0.584-01 


F2 


-0.369-00 


-0.340-00 


-0.165-00 


-0.665-02 


N.R. 










Fi 


-0.032-01 


-0.296-01 


-0.151-00 


-0.550-01 


F2 


-0.369-00 


-0.342-00 


-0.174-00 


-0.636-02 


P.F. 










Fi 


-0.101-01 


-0.372-01 


-0.140-00 


-0.283-01 


F2 


-0.324-00 


-0.293-00 


-0.119-00 


-0.022-02 



Results for an inelastic transition are given in Table 2. The results obtained in the 
non-relativistic approach compare well with the exact ones over the full range of Q^. 
The slight discrepancies are quite similar to those observed for the elastic form factors. 
While the point-form results compare well with the other ones in the intermediate range, 
0.1 < Q'^ /m? < 3, they fail at low and at large Q^. In the former case, Fi{q^) does not 
go to zero when 0, violating the current conservation, Eq. @). In the latter case. 
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F^io^) evidences a change in sign around ~ 10 m^, also preventing one from fulfilling 
this relation. This can be traced back to the different behaviour of the intermediate form 
factors Fi and F2, which keep the same sign but scale like and at high Q^, 
respectively. These results for an inelastic transition complement those for the elastic 
case. Again, a large discrepancy with exact results appears, but the bad behaviour of 
form factors at low as well as large momentum transfers is more clearly correlated with a 
violation of current conservation. 



4 Discussion and conclusion 



The present study was motivated by the necessity to check the reliability of retaining the 
single-particle current for the calculation of form factors in the point form of quantum 
relativistic mechanics, which was recently employed in different works p, |^. With this 
aim, we considered a simple model where the exact result is known. It was found that the 
non-relativistic approach does particularly well up to momentum transfers of 3-4 times the 
constituent mass, including the case of a strong binding. While the point-form approach 
does correctly for small bindings and small couplings, large discrepancies with the exact 
results appear as soon as the momentum transfer or the coupling increases. Detailed 
examination evidences three features. 

- The form factors decrease more rapidly than they should. This points to an extra ^ 
dependence that is absent in the exact calculation. 

- For strong couplings, corresponding to a sizeable presence of high momentum compo- 
nents in the wave function, the charge radius turns out to be much larger than the exact 
one. 

- Finally, as emphasised by the results for an inelastic transition, current conservation is 
strongly violated. 

Although it is not quite certain, there is good reason to believe that the non-relativistic 
calculation does relatively well because it fulfils current conservation. With this respect, 
the failure of the point-form approach is likely to reside in the incomplete character of 
the current operator. Current conservation may be enforced by the replacement 

(jM) ^ (JO - W ■ gV?'- 

Apart from the unsatisfactory character of this recipe, due to the presence of a pole at 
= , it does not solve the problem of the too fast drop-off of the elastic charge form 
factor at high Q^. Most probably, two-body currents have to be considered. This is 
consistent with the known fact that, in the point form approach, the interaction appears 
in all the components of the 4-momentum operator |jl5[. Whether a minimal set of two- 



body currents will be sufficient to provide results in better agreement with the exact ones 
is not clear however. These currents should also correct for the failure of the point-form 
results to reproduce the Born amplitude expected from the underlying field-theory model. 

The above results cannot be applied directly to the calculation of the deuteron or 
nucleon form factors @, 0. However, the qualitative similarity in the results strongly 
suggests that the kinematical boost, which provided a nice description of the nucleon 
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form factors, represents an incomplete account of relativistic effects. More likely, the 
agreement is the consequence of neglecting significant contributions to the current. This 
conclusion is to be preferred with two respects. 

- It leaves some room for the well known contribution of the coupling of the nucleon to 
the photon through vector-meson exchange (vector-meson dominance mechanism) , which 
roughly provides half of the proton's squared charge radius. 

- On the other hand, it leaves room for another relativistic effect related to the nature of 
the coupling of the constituents to the exchanged boson. This effect increases the form 
factor at high rather than the opposite as in ref. [^. It is known to explain the different 
asymptotic form factors in a non-relativistic and a relativistic calculation in QCD, which 
scale like and respectively (see for instance refs. |TB|, P^). 



Altogether, two-body currents should produce quite sizeable contributions. Their role 
seems to be more essential in the point form than in other approaches. 
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